Quantum magnetism and counterflow supersolidity of up-down bosonic dipoles 
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We study a gas of dipolar Bosons confined in a two-dimensional optical lattice. Dipoles are 
considered to point freely in both up and down directions perpendicular to the lattice plane. This 
results in a nearest neighbor repulsive (attractive) interaction for aligned (anti-aligned) dipoles. 
We find regions of parameters where the ground state of the system exhibits insulating phases 
with ferromagnetic or anti-ferromagnetic ordering, as well as with rational values of the average 
magnetization. Evidence for the existence of a novel counterflow supersolid quantum phase is also 
presented. 
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I. INTRODUCTION 

Recently, the physics of ultracold dipolar gases has re- 
ceived growing attention. Experimental and theoretical 
studies are motivated by the long-range and anisotropic 
dipole-dipole interaction which introduce a rich variety 
of quantum phases For instance, supersolid and 

checkerboard phases are predicted in the phase diagram 
of polarized dipolar Bosons confined in a two-dimensional 
(2D) optical lattice Q. In bilayer samples where dipole 
interactions can also be attractive more exotic quantum 
phases are accessible, for example a pair-super solid [3[. 

Experiments with ultracold atoms have now high- 
lighted the role of dipole-dipole interactions. Notably, 
with 52 Cr, control via Feshbach resonances allows one to 
efficiently reduce contact interactions and enter a regime 
where magnetic dipole-dipole interactions become domi- 
nant [4]. Polar molecules also appear very promising to- 
wards implementation of degenerate dipolar gases (see 
and References therein), in particular since the demon- 
stration of ultra-cold rubidium-potassium, and lithium- 
cesium molecules prepared in their ground rotovibra- 
tional state [H, @] . 

Note that in most cases, if not all cases considered so 
far, dipolar gases are polarized, i.e. magnetic or elec- 
tric dipoles point in the same direction. This, however, 
does not have to be the case always. A prominent exam- 
ple is that of exciton gases. Excitons, which are bound 
electron- hole pairs, obviously carry an electric dipole mo- 
ment that can a priori attain quite arbitrary directions. 
Nevertheless, in indirect quantum wells, electrons and 
holes are confined in spatially separated regions such that 
electric dipoles can be aligned [7j or even anti-aligned in 
type-II heterostructures [8|. 

Another extreme case concerns molecules that follow 
the Hund's rule (a) [9j]. For such molecules the electric 
dipole is either parallel or anti-parallel to the direction 
of the magnetic moment. If one takes thus a sample of 
such molecules and polarizes it magnetically, say in the 
up direction, one will obtain in this way a dipolar gas 
with certain fraction of electric dipoles pointing up, and 
the remaining fraction pointing down. This is the situa- 
tion that we study in the present paper. In fact there was 
recently a spectacular progress in cooling and trapping 



of magnetically confined neutral OH molecules. This 
progress opens important perspectives to study novel 
quantum phases with ultracold dipoles [To|, [ll|. In its 
ro-vibronic ground state OH is indeed a Hund's case (a) 
molecule, and it features both an electric d = 1.67 D and 
a magnetic dipole moment \i — 1.2/i B (with /i B being the 
Bohr magneton), which can be aligned or anti-aligned 
independently. Using these two degrees of freedom, OH 
molecules may be confined in a 2D optical lattice, where 
the dipolar interactions lead to rich quantum phases, as 
we show in this work. 

We study here a sample of bosonic dipoles confined 
in a 2D lattice. We consider that dipoles are free to 
point in both directions normal to the lattice plane. This 
constitutes a novel ingredient compared to our previous 
works [III, [l3[ , and results in a nearest neighbor interac- 
tion either repulsive for aligned dipoles, or attractive for 
anti-aligned ones. We consider the case of dipolar inter- 
actions to be relatively weak compared to the contact in- 
teractions, and remarkably we find that the dipole-dipole 
interactions dominate the physics of the system. Using 
a mean-field Gutzwiller approach, we show that the sys- 
tem presents Mott-insulating phases with ferromagnetic 
or anti-ferromagnetic ordering, as well as with fractional 
values of the average magnetization, depending on the 
imbalance between the population of dipoles pointing in 
opposite directions. We also find evidences of a novel 
counterflow supersolid quantum phase. The latter ex- 
hibits broken translational symmetry, namely a modula- 
tion of the order parameter on a length scale larger than 
the one of the lattice spacing, analogously to supersolid 
phases. 



II. HAMILTONIAN OF THE SYSTEM 



We consider a sample of dipoles in the presence of a 2D 
optical lattice, and an extra confinement in the perpen- 
dicular direction. The dipoles are free to point in both 
directions normal to the lattice plane, as shown in Fig. [TJ 
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p 2 (k) are the Fourier transform of the dipole potential 
and the density respectively. We consider that the s- 
wave scattering length is independent on the orientation 
of the dipoles. Instead, the on-site dipolar contribution 
Udd depends both on the orientation of the dipoles and 
on the geometry of the trapping potential, and it can 
be varied by changing the ratio between the vertical to 
the axial confinement. For simplicity we will focus on 
the specific case of a spherically symmetric confinement, 
where the on-site dipolar interactions average out to 
zero Udd = 0, and the resulting on-site interactions are 
all equal to U. We consider the case of dipole-dipole 
interactions to be 600 times weaker with respect to the 
on-site interaction, i.e. /7nn = U/600 [lij . 



FIG. 1: Schematic representation of a 2D optical lattice pop- 
ulated with dipolar Bosons polarized in both directions per- 
pendicular to the lattice plane. The particles feel repulsive 
intra-species U aa , Ubb, and inter-species U a b repulsive on-site 
energies. The first nearest-neighbor interaction is repulsive 
£/nn > for aligned dipoles, while it is attractive — £/nn for 
anti-aligned particles while the hopping term J is equal for 
both the species. 



The system is described by the Hamiltonian 
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where a = (a, b) indicates the type of species, i.e. dipoles 
pointing in the up and down direction perpendicular to 
the 2D plane of the lattice, respectively. U aa an d Ubb 
are the on-site energies for particles of the same species, 
while U a b is the on-site energy for different species. The 
long-range dipolar interaction potential decays as the in- 
verse cubic power of the relative distance r^, which we 
express it in units of the lattice spacing. For computa- 
tional reasons, in most theoretical approaches the range 
is cutoff at a certain neighbors. In the present work, 
we consider a range of interaction up to the 4th nearest 
neighbor. The first nearest neighbor dipolar interaction 
is attractive (repulsive) for particles of the same (dif- 
ferent) species, with strength Unn > (or respectively 
— Unn), and we consider an equal tunneling coefficient 
(J) for both the species, while the densities of the dipoles 
pointing upwards and downwards are fixed by the corre- 
sponding chemical potentials fi a . 

The on-site interactions have two contributions: 
one is arising from the s-wave scattering length 
U 8 = 47r ^ Qs / d 3 rp(r) 2 , and the second one 
is due to the on-site dipole-dipole interaction 
U ^ = (2^F f d3kU dd(k) p 2 (k), where U dd (k) and 



A. Filling factor and imbalance 

The properties of the system are conveniently ex- 
tracted using the operators given by the sum (filling fac- 
tor) and by the difference (imbalance) of the two species 
number operators at each site of the lattice, namely by 
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which are simultaneously diagonal on a given Fock state 
\v,m)i. Notice that the eigenvalues of these two opera- 
tors are not independent. In fact, by fixing v the eigen- 
values of rhi can only assume 2^ + 1 values given by 
m = {— z/, — v + 1, -f^}, in complete analogy with the 
angular momentum operator Sf and its projection along 
the z axis Sf, as we will discuss in Sec. Ill Bl It is also 
useful to introduce the average magnetization of the sys- 
tem, defined as 
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where N s is the total number of lattice sites. 

Substituting Eqs. (|2]) into equation (pQ) allows us to 
express the system Hamiltonian as H = Hq + H™ + H^ m ? 
where 
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In Eqs. dHHJ), we have introduced the chemical potentials 

Ma i l^b 

M± = ^ ' ^ 

which respectively fix the eigenvalues of the filling factor 
(+), and the imbalance operators (— ) in Eq. ©. In the 



following we consider H^ m to be a small perturbation 
on the interaction terms ((4]) and ((5]). In the limit where 
U ^> (Unn,J), the ground state of the system is found to 
be a uniform distribution of constant filling factor V{ = v 
at each site of the lattice. The value of v is fixed by 
and can be integer as well as semi-integer. This is 
better understood at J = 0, where we can calculate the 
expectation value of Hq on a given classical distribution 
of atoms in the lattice |3>) = Yii Wi,rrii)i, as follows 



(8) 



where V{ = (3>|z^|<I>). In the right hand-side of Eq. (|8j) 
all sites i are equal, and like in the homogeneous case 
of a Bose-Hubbard Hamiltonian at J = [15], the min- 
imum of Eq. (|Sj) is provided by a uniform distribution 
in the lattice, where V{ = v at each site. Instead, for a 
given z/, finding the distribution of = ($|m^|$), which 
minimize the expectation value 
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FIG. 2: Ground state of the system at J — 0, calculated for a 
4x4 elementary cell satisfying periodic boundary conditions, 
and ?7nn = £7/600. The text in parentheses (y,M) indicates 
the filling factor v and the average magnetization M. 
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is non-trivial due to the presence of the long-range in- 
teraction. However, we can qualitatively argue that for 
^> Unn, the minimum of the energy ((9]) is ob- 
tained for mi = v x sign[/i_] Vi, which corresponds to 
a ferromagnetic phase (FM) of average magnetization 
M = v x sign where only particles of one species 
are present. Instead, for fi_ = 0, a succession of near- 
est neighbors with ra$ = v and rrij = —v provides the 
minimum of Eq. (j9j), and the phase is anti- ferromagnetic 
(AM), i.e. M — 0. In other words, the spatial distribu- 
tion of the particles is given by sites occupied from the 
species a alternated with sites occupied by the species b 
in a checkerboard-like structure. In Fig. [2] we plot the 
ground state phase diagram at J = 0, in the \i_ vs. /i + 
plane, where the text in parentheses indicates filling fac- 
tor and average magnetization, i.e. v and M respectively. 
Note that between the AM, and the FM phases, we find 
magnetic phases which present non trivial rational values 
of the average magnetization (RM). The precise choice of 
the cutoff range in the dipolar interaction, and the lat- 
tice size, determine the fractional character of the allowed 
ground state magnetization. 



In the next section, we include the presence of tunneling. 
The physical properties of the system are described by 
an effective Hamiltonian, supported by the existence of a 
low-energy subspace. Hence, the theoretical description 
of the system cannot be based on standard mean-field 
theory, which is not suitable to correctly describe the 
ground state of the system, as we discuss in the following 
section. 



B. Low-energy subspace and effective Hamiltonian 

The ground state of the system at J = is described 
by a product over single-site Fock states of the type 
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with uniform total on-site occupation 2v. As single par- 
ticle hopping changes the total on-site population, it 
breaks the translational invariance of the ground state 
with respect to the total on-site occupation 2v. The en- 
ergy cost of these excitations is of the order of the on-site 
interaction energy J7, and is therefore very costly in the 
limit where U ^> (Unn,J)- On the contrary, exchang- 
ing two particles from nearest neighboring sites does not 
require such a large amount of energy. This defines a 
low-energy subspace spanned by the \a) configurations at 
constant filling factor v of Eq. ([TO]) , which is energetically 
well separated from the rest of the Hilbert space in the 
limit of parameters we consider. Thus, a successful de- 
scription of such a system is obtained through an effective 
Hamiltonian H e ^ restricted to the low-energy subspace, 
where single-particle hopping is suppressed and tunnel- 
ing is included at second order in perturbation theory. 
The validity of the effective Hamiltonian relies on the 
existence of this low-energy subspace well separated in 
energy from the subspace of virtual excitations, to which 
it is coupled via single-particle hopping. The relevant 
virtual subspace is then obtained from the states \a) via 
single particle hopping, and it is spanned by the states 
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as schematically represented in Fig. El 
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For a given the second line of the Hamiltonian (fT3|) 
can be equiyalently written in terms of the spin operators 
at site i 

Ham, 

given by 



(16) 
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FIG. 3: Schematic representation of a two-particle hopping 
between the states \ot) and belonging to the low-energy 
subspace at v — 1/2. These states are coupled through virtual 
excitations to the states I7) by single-particle hopping. 



This situation is in fact similar to the one discussed in [3| 
for a bilayer optical lattice, and therefore we apply the 
same technique to compute H e ^ . In the basis of constant 
total on-site population 2^, the matrix elements of such 
a Hamiltonian at second order in perturbation theory are 
given by 

{a\H eS \/3) = (a\H \P)-±Y,( a \fir\l)(l\Hr\a) 



(12) 
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where Ho = Hq+H™, given by the sum of the interaction 
terms (jlj and ([5]), is diagonal on the states |a), and the 
single-particle tunneling term H" m of Eq. (j6]) is treated 
at second order. For a given state \a), 



E a = U + U NN Am^ N , 



(13) 



with AmJJ N 



E/c^i 2rn k /\r ik \ 3 - J2 k ^j 2rn k /\r jk \ 3 - 1, 
where rrii indicates the population imbalance at site i of 
Eq. (|2j). For U ^> Unn, the denominators E 1 — E a are 



all of order J7, which leads to 
2J 2 v- 
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where q = a^J" and = a\bi are composite operators, 
corresponding to the creation of a particle of one species 
and a hole of the other species, such that 



Ci\v,mi) \/v(y + 1) - miimi - \)\v, m, - 1) 
c\\v,rrii) = \/v{y + 1) - m,(m, + l)|z/, m, + 1), 

while their commutation relation obeys [Q,ct] 
—2rhiSij. 



(15) 



where 3 = {cF Xl a yi a z ) are the Pauli matrices, and 
u = (a, b) indicates each species. Thus, the creation and 

annihilation operators (fT5j) become c\ = Sf + zS'f and 
q = 5f — z^f respectively, while the imbalance operator 
is given by rhi = Sf as already anticipated in Sec. Ill Al 
Therefore, in a spin representation as in Eq. ([T6]h the 
second line of the Hamiltonian (fH|) acts as a Heisenberg 
spin Hamiltonian (see e.g. [16]). The chemical potential 
/i_ plays the role of an external magnetic field along the 
z axis. The interplay between and the long-range 
interaction determines the magnetic ordering of the sys- 
tem, as we will discuss in the next section. 



III. MEAN-FIELD 

In this section, we provide a mean-field solution to the 
effective Hamiltonian (fT4|) to investigate quantum phases 
of the system. We identify the different phases through 
the composite order parameters (cj), as well as both the 
single-particle ones (a.j), and (bi). 

For every subspace at constant filling, we find that 
the system presents three different kinds of phases. The 
Mott-insulating phase (MI), with a well defined number 
of particles at each site of the lattice, and absence of any 
low-energy transport [15]. The MI is characterized by 
vanishing (q) = (a$) = (bi) =0, and depending on the 
value of fi_ presents FM, AM or RM ordering. The sec- 
ond quantum phase, a super- counter- fluid phase (SCF), 
exhibits on-site density fluctuations where a net trans- 
port of atoms is suppressed but a counterflow is present, 
for which the currents of the two species have equal am- 
plitudes but opposite directions [17j. In the SCF phase, 
while the single-particle order parameters still vanish 
(hi) = (bi) =0, the composite order parameters are non- 
zero (di) 0, indicating the presence of counterflow. We 
also find evidence of a third and novel quantum phase, 
namely a counterflow supers olid phase (CSS). The CSS is 
characterized by vanishing single-particle order parame- 
ters (hi) = (bi) =0, and non- vanishing composite order 
parameters (q) 7^ 0, coexisting with broken translational 
symmetry, namely, a modulation of both m^, and (q) on 
a scale larger than the one of the lattice spacing, analo- 
gously to the supersolid phase. 

To determine the insulating phases we perform a per- 
turbative treatment at first order in the composite or- 
der parameters ipi = (cj), which allows us to compute 
the boundaries of the insulating lobes. Furthermore, we 
solve the time dependent Gutzwiller equations in imagi- 
nary time to determine the nature of the SCF-CSS phases 
outside the lobes. 
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A. Insulating lobes 

The low-energy subspace is spanned by the classical 
distribution of atoms in the lattice \a) of Eq. ([TO]) . Simi- 
larly to the two layer system discussed in [3| , in the limit 
where U ^> C/nn? asymptotically all classical states \a) 
become stable with respect to single-particle-hole exci- 
tations and develop an insulating lobe at finite J. The 
energy of single particle-hole excitations is of the order 
of U at J = and is given by the width of the lobes at 
finite J (see, e.g., the thin blue/black lobes in Fig. [4]). 

Instead, the low-lying excitations remain within the 
subspace and are obtained by adding (PH) or removing 
(HP) one composite, made of a particle of the upper- 
polarized dipoles (species a), and a hole of the lower- 
polarized dipoles (species 6), at the i-th site of the lattice. 
This corresponds to flipping the direction of a dipole at 
the site z, respectively from down to up (PH) or from 
up to down (HP). For any given configuration |a), one 
can calculate the corresponding energy costs using the 
diagonal terms of the effective Hamiltonian (fl4]h which 
are respectively given by 

EUJ) = -2 M _ + 4Z7 NN Tpti ~ Hr E m ^ 



(17) 
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Note that in the last expressions, there is no explicit de- 
pendence on the chemical potential This is because 
by adding or removing one composite, we remain within 
the subspace at filling factor ^, and therefore the contri- 
bution of fi + vanishes in the calculation of the excitations 
([TT]) . By using a perturbative mean- field method, we can 
calculate the order parameters ipi = (q) for |a), which 
satisfy the equations 
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where ipi = ^2/j\. ^j- With Eqs. ([T8]) one can calculate 
the mean-field lobes of any distribution of atoms in the 
lattice |a), provided that the elementary excitations (fT7|) . 
are positive in some parameters range. For every site i 
of the lattice, one can adequately flip the direction of a 
dipole, and depending on mi either one of the following 
conditions apply 



/i_ < 2C/NN 1 T7 

^ hj\ 3 



k^i 



/i_ > 2/7 N n^ 



rrij 



(19) 



For example, suppose the site i is occupied only by one 
particle of the species a, i.e rrii = 1/2, then for this site 
the conditions (fT9j) reduce to /i_ > 2C/nn Sj^i m j'/l r ij| 3 5 



since the only possible excitation at this site corresponds 
to remove a composite. Conditions ([T9]) are necessary 
for the existence of an insulating lobe and provide its 
boundaries at J = 0. For each site of the lattice one 
has such a condition (Qj)j), and an Eq. ([T8]) . The latter 
constitutes a set of coupled equations, which can be writ- 
ten in the matrix form A^(/i_, U, [/nn, J) ■ ^ = 0, with 
x/j = (• • • ipi • • • ) T being the vector of the order param- 
eters at each site of the lattice. For every a non- 
trivial solution is provided by the smallest J for which 
det[7W(/i_, U, £7nn, J)] — 0, that is the insulating lobe of 
the \a) configuration in the J vs. \i_ plane. In Fig. [5] 
we plot the ground state insulating lobes calculated in 
this way for v = 1/2 (left) and v = 1 (right). For all 
filling factors z/, we find an AM ground state (z/, M = 0), 
which presents a spatial distribution of alternating sites 
occupied by particles of species a and 6 resembling a 
checkerboard structure. Remarkably, the larger the ^, 
and the more stable is the AM ordering with respect to 
flipping the direction of a dipole. Indeed, using the in- 
equalities ([T9]h it is not difficult to calculate the bound- 
aries at J = of such a checkerboard structure, which, 
for a fixed v are given by — IzvtUm* < fJ>~ < 2zv£Unn, 
where £ = 1 - 2~ 3 / 2 - 2~ 3 + 2 x 5~ 3 / 2 , and z = 1 
is the coordination number (here z = 4). Instead the tip 
of the AM lob es is found to be ^-independent, and given 
by J/U = \J £ £/nn/2£7 at /i_ = 0. By increasing the 
absolute value of fi_ we find RM ground states with ra- 
tional values of the average magnetization, corresponding 
to M = (zL2i/, ±4^, d=6i/)/8. The exact fractional values 
of M in the ground state, depend on the cutoff range 
in the dipolar interactions, and on the size of the lat- 
tice. We have used a 4 x 4 elementary cell with periodic 
boundary conditions, and dipolar interaction range cut 
at the 4th nearest neighbor. By considering more neigh- 
bors in the interactions, and larger lattices, we expect to 
find RM states appearing at all rational M, asymptoti- 
call y ap proaching a Devil's staircase as recently shown in 
pH.l22f. Finally we find a FM ground state (y, M = ±v), 
in which only particles of one type are present. 

It is worth noticing that the insulating lobes calculated 
in this way, do not contain any dependence on which 
does not enter into Eqs. ([T8]) as previously discussed. 
Therefore, for any given value of in order to obtain 
the ground state phase diagram one has to compare the 
energies of the ground state configurations at different 
v. Using the effective Hamiltonian ([T4]h for any value of 
J, and we compare the energies of the ground 
state configurations for different z/, and select the state 
with the smaller energy. In this way we have obtained 
the phase diagram at J = of Fig. [21 



B. Counterflow superfluid/supersolid 

In the low-energy subspace at constant z/, the 
Gutzwiller Ansatz on the wave function of the system 
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FIG. 4: Insulating lobes at v — 1/2, for fi + — OAU, and 
£/nn = £7/600. The thick red line is the ant i- ferromagnetic in- 
sulating state calculated with the effective Hamiltonian . 
The thin blue lines and the black dashed lines represent in- 
sulating lobes calculated with respect to single-particle-hole 
excitations of the species a and species b respectively, and are 
some of the many dominant configurations of the J = 0.035?7 
and ii- = point (cross). The two species are sketched as a 
plain dots (a) and squares (b). 



reads 



(20) 



i m= — v 



where we allow the Gutzwiller amplitudes to depend on 
time, i.e. fu%.(t). We obtain the equations of motion for 
the amplitudes by minimizing the action of the system, 
given by S = J dtC, with respect to the variational pa- 
rameters fv)m(t) and their complex conjugates /^,m(£), 
where 



is the Lagrangian of the system in the quantum state |<E>) 
(23j . Therefore setting to zero the variation of the action 
with respect to f^m leads to the equations 

4J 2 



ih— f (i) 



2/i_ 



U 



(oh 



I 13 " L Jv,m 



+ ^v^+i)-^+i)/l;i4 ( 22 ) 

where (m,> = Em=-^ m l/^m| 2 , the fields ^ = i/jj, 
J2(j}.(rhj), and Ylj^ii^j) /\ r ij\ 3 nave to ^ e calculated 



in a self consistent way, and the order parameter ipi = 
($|£i|3>) is given by 

V 

^ = E >(" + l)-m(m + l)/;«/;i +1 . (23) 



We solve Eqs. ([22]) in imaginary time r = it, which due to 
dissipation is supposed to converge to the ground state. 
In Figure [5] we show the ground state phase diagram of 
the system for v = 1/2 (left) and v = 1 (right), com- 
puted in this way for 17 NN = U/600. The phase diagram 
is symmetric with respect to the fi_ = axis. This is 
because the description of the system is identical under 
the interchange of the species a with the species 6, and 
vice versa. For v = 1/2, in the region immediately out- 
side the insulating AM lobe, depending on the values of 
J and \i_ we find either SCF or CSS. Note that we don't 
find any evidence ofCSSat/i_=0, which indicates that 
a finite imbalance between the two species is a necessary 
condition for the system in order to sustain CSS. 

Let us underline that for v = 1/2, our effective Hamil- 
tonian ([T4]) can be mapped onto a hard-core dipolar Bose- 
Hubbard Hamiltonian, provided that we neglect a con- 
stant factor proportional to and rescale the first near- 
est neighbor interaction by the small quantity —2J 2 /U. 
Contrary to mean-field predictions, Monte Carlo studies 
of the hard-core Bose-Hubbard Hamiltonian on square 
lattices, with interactions extended to first and second 
nearest neighbors, have shown that no super solid phase is 
obtained by doping the checkerboard solid (AM ordering) 
[l9|,[20|. However, it was recently demonstrated by Monte 
Carlo simulations that considering an infinite range in 
the dipolar interactions stabilizes the supersolid phase, 
obtained by doping the checkerboard solid (AM) either 
with particles or holes [21], in agreement with mean- field 
predictions. In our treatment we are approaching this 
limit, and we therefore expect the CSS phase to be sta- 
bilized by the long-range interactions. Furthermore, be- 
cause J > 0, the rescaling of the first nearest neighbor 
interaction by the small quantity —2J 2 /U in the map- 
ping process, should enforce further this effect. Although 
our treatment is approaching infinity in the interaction 
range, it would be important to verify our predictions 
with first-principle quantum Monte Carlo simulations. 

To get reliable results, one should combine the Gutzwiller 
predictions with an estimate of the limits of validity of 
i?eff ' beyond which the subspace of constant v looses its 
meaning. Before starting the discussion on the validity of 
the subspace, let us explain how we define the dominant 
classical configurations of a given state |3>). It is not 
difficult to see that Eq. ([20]) can be equivalently written 
as 



(24) 



{rh} 



where m = (mi, m^, m N ) is a collection of the in- 
dices m at each site, and we have introduced the notation 
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FIG. 5: Ground state of a 4 x 4 square lattice satisfying pe- 
riodic boundary conditions, for v = 1/2 at fi + = 0.5U (left), 
and v = 1 at /x+ = 1AU (right), and Unn = U/600. The 
text in parentheses (zy, M) indicate the filling factor is, and 
the average magnetization M, respectively. 

9 fa — Eli /^mi • The advantage of writing the Gutzwiller 
state |<E>) in the form lays in the product over single- 
site Fock states \a) — Y\i\ v i m i)ii which is nothing but 
a classical distribution of atoms in the lattice. Therefore 
we can rewrite Eq. ([24]) as 

|*> = £> a |a>. (25) 

{a} 

For each point of the phase diagram, from the ground 
state Gutzwiller wavefunction, we define the domi- 
nant classical configurations with the criteria \g^,\ — 
irii/^mj > (0.02) 4 , and we require l/^il 2 > 0.05, im- 
plying that each of the contributing fj>%n should also be 
sufficiently large. For each of these configurations, we 
calculate the lobe with respect to single particle-hole ex- 
citations [27]. If the system at this given point of the 
phase diagram turns out to be stable against all dom- 
inant single particle- hole excitations (in other words, if 
this point is inside all selected single particle- hole lobes), 

H^fi is considered valid. This procedure is shown for 
/i + = 0AU, J = 0.035/7 and ji_ = in Fig. U and gives 
the thick black lines of Fig. (j5j). On the right hand side 
of these black lines, the low energy subspace is not well 
defined and therefore the effective Hamiltonian looses its 
meaning, leaving the description of the system to the do- 
main of single-particle single-hole excitation theory that 
predicts SF and SS phases for each component separately. 

We have already mentioned, that the boundaries of the 
lobes calculated with the effective Hamiltonian do not 
show any dependence on the chemical potential ja + , which 
does not give any contribution in the expression of the 




FIG. 6: Slices of the 3D ground state of a 4 x 4 square lat- 
tice satisfying periodic boundary conditions, for v — 1/2 and 
C/ NN [7/600. 



low- lying excitations ([17]) . This is not true in the case 
of the single-particle-hole insulating lobes, since adding 
or removing a single particle results in a change of both 
and This makes the process of estimating the 
limits of validity of more complicated and leads to 
a 3D phase diagram in J, and \i_ highly non-trivial. 
In Fig. [6] we present slices of the 3D phase diagram, 
calculated at constant values of for v = 1/2 only. 
From Fig. [6j it is evident that the insulating magnetic 
phases persist for a wide range of ja + values. Instead, 
the SCF and the CSS phases are more affected by the 
limits of validity of , which may imply high degrees of 
control in order to experimentally observe this quantum 
phases. 

Finally, it is useful to estimate the critical tempera- 
tures at which we expect these quantum phases to be ex- 
perimentally observable. Bose-Einstein condensation in 
alkali atoms occurs at temperatures of the order ~ 100 
nK. Assuming that condensation of OH molecules oc- 
curs at similar temperatures we expect the SCF to be 
observed at ~ 100 nK. For the AM insulating phases, 
the gap at J = is of the order ~ 2v x 10 -2 Z7, and 
the critical temperature is expected to be of the order 
of the standard insulator-superfluid transition tempera- 
ture, i.e. 10 to 50 nK [24]. Therefore, the CSS phase is 
expected to be observed somewhere in between 50 to 100 
nK. The lowest measured temperature of ultracold atoms 
in optical lattices is ~ 1 nK [25[ , where superexchange in- 
teractions, of the order ~ J 2 /U, are dominant. Recently, 
a two-component hard-core Bose-Hubbard Hamiltonian 
was investigated with Monte Carlo techniques [26], and 
insulating AM as well as SCF phases where observed. In 
[26[ , the largest critical temperature observed for the AM 
phase is T c - 0.12J, at J/U = 0.0125, and T c - 0.104J 
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at U/J = 13 for the SCF phase. 

IV. CONCLUSIONS 

In conclusion, in this work we have studied a gas of 
dipolar Bosons confined in a two-dimensional optical lat- 
tice, where the dipoles are free to point in both directions 
perpendicular to the lattice plane. We have found regions 
of parameters where the ground state presents insulating 
phases with antiferromagnetic or ferromagnetic ordering, 
as well as with rational values of the average magnetiza- 
tion. Our mean-field calculations predict the existence 
of a novel counterflow super solid quantum phase, which 
presents a spatial modulation of the density coexisting 
with the presence of counterflow. Our work in a sense 
is the first step in the studies of dipolar gases with non- 



polarized dipoles, and is relevant for experimental studies 
of the ultracold OH molecules. We expect that the meth- 
ods and ideas developed here will turn out to be useful 
also for the study of more complex systems, such as ex- 
citons in indirect quantum well structures, or completely 
unpolarized dipolar gases. 
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